Dumbbell diffusion in a spatially periodic potential 
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We present a numerical investigation of the Brownian motion and diffusion of a dumbbell in a 
two-dimensional periodic potential. Its dynamics is described by a Langevin model including the 
hydrodynamic interaction. With increasing values of the amplitude of the potential we find along 
the modulated spatial directions a reduction of the diffusion constant and of the impact of the 
hydrodynamic interaction. For modulation amplitudes of the potential in the range of the thermal 
energy the dumbbell diffusion exhibits a pronounced local maximum at a wavelength of about 3/2 
of the dumbbell extension. This is especially emphasized for stiff springs connecting the two beads. 

PACS numbers: 87.15.Vv, 82.35.Lr, 05.40.-a 



Introduction.- Investigations on the diffusion of dif- 
ferent colloidal particles in a homogeneous solvent have 
a long history [l|, |2|, while studies on the diffusion of 
small spheres, dimers and polymers in different poten- 
tials attract considerable interest only for a short time 
[1) [1, IE 0) 0, B Hi • Laser- tweezer arrays are a new pow- 
erful tool for generating the desired spatially periodic, 
correlated or unstructured potentials in order to study 
the effects of inhomogeneous potential landscapes on the 
motion of colloidal particles [3, 0, B EB] • Furthermore 
recent studies of dumbbells and polymers in randorn po- 
tentials are exciting issues in statistical physics [III Sjl ■ 

Several of these investigations are motivated by possi- 
ble applications like particle sorting in inhomogeneous 
potentials. For example cross-streamline migration of 
colloidal particles has been found in a flow through an 
optically induced periodic potential. Since this migra- 
tion depends on the extent of the colloidal particles, the 
laser- tweezer array has been successfully used for sorting 
particles with respect to their size [1, [To] . 

We investigate the Brownian motion and the diffusion 
of dumbbells through a two-dimensional periodic poten- 
tial, which is described by a Langevin model. In doing 
so we include the effects of the hydrodynamic interac- 
tion between the two beads of the dumbbell and focus 
on the interplay between the dumbbell extension b and 
the wavelength A of the spatially periodic potential. In 
the context of this work the dumbbell may be considered 
as a simple model for several deformable objects such as 
pom-pom polymers fl3| or two colloids which are con- 
nected either by a rather flexible A-DNA molecule or by 
a semi-flexible actin filament. 

Our numerical studies reveal a significant dependence 
of the dumbbell diffusion on the ratio A/6. With the 
potential amplitude Vb of the order of the thermal ex- 
citation energy fc^T we find a remarkable maximum of 
the dumbbell diffusion constant in the range of A « 36/2, 
whereby the height of the diffusion maximum increases 
with the stiffness of the spring connecting the two beads 
of the dumbbell. Another remarkable effect is the re- 
duction of the influence of the hydrodynamic interaction 
with increasing potential amplitude. 



Model.- We describe the Brownian motion of a dumb- 
bell in a two-dimensional periodic potential by the fol- 
lowing Langevin equation (without inertia) 
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for the bead positions r; = {xi,yi, Zi). The linear spring 
force i^j between them is determined by the harmonic 
potential , 



$(ri,r2) 



(6-|ri-r2|)^ 



(2) 



with the equilibrium length b of the spring and the corre- 
sponding spring constant k. The spatially periodic force 
FY = — Vy(ri) is derived from the two-dimensional po- 
tential in the xy-planc 
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which can be realized in experiments by a laser tweezer 
array. Its amplitude Vq may be changed by varying the 
intensity of the laser light. The same wavelength A is 
chosen in the x- and in the y-direction. 

In the absence of hydrodynamic interaction (HI) be- 
tween the beads the mobility matrix H is a diagonal ma- 
trix (Hii = ^ I , Hij = for i ^ j) being inversely propor- 
tional to the Stokes friction coefficient ( — Gnrja which 
depends on the solvent viscosity rj as well as on the ef- 
fective hydrodynamic bead radius a. The HI between 
the two beads is taken into account by the Rotne-Prager 
tensor [l4[ where the mobility matrix for i ^ j has the 
following structure: 
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Note that 



To is the distance vector between the 



beads and is its norm. 

The stochastic forces Ff caused by the thermal heat 
bath are related to the dissipative drag by the fluctua- 
tion dissipation theorem which ensures the correct equi- 
librium properties. They can be combined to a single 
supervector F'^ = (Ff 



, F2 ) that reads 



(5) 
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T is the temperature, ks the Boltzmann constant and 
^{t) is the uncorrelated Gaussian white noise vector with 
zero mean and unit variance 

m) = 0, (6) 
{me{t')) = 5{t~t')\. (7) 

The fixed parameters in our simulations are 6=1 for 
the typical length scale and fc^T = 1 which determines 
the energy scale. If not stated otherwise we use f = I, 
rj = \. k determines the binding energy of the spring in 
units of the thermal energy . Most of our results were 
observed after averaging over more than 10* ensembles. 

Results.- A typical trajectory of the center of mass of 
the dumbbell (CM), R = (ri + r2)/2, in the a:y-plane is 
shown in Fig. □ for fc = 10, A = 2& and Vb = 2fcBT. It 
passes the saddles between the maxima of the potential 
and accordingly the trajectory adapts to the quadratic 
structure of the potential landscape. If one of the two 
parameters k or \/h is reduced, larger excursions of the 
CM away from the potential minima occur and even di- 
agonal jumps between the valleys are found frequently. 
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FIG. 1: The trajectory of the CM of the dumbbell in the xy- 
plane follows predominantly along the saddles between the 
minima of the potential. Parameters: k = 10, Vo/kBT — 2, 
and A = 26. 

The mean square displacement {Rf{t)) — 2Dit , {I — 
X, y, z) increases linearly in time along each spatial di- 
rection as shown for one parameter set in Fig. [2l This 
behavior is typical for normal diffusion. Parallel to the 
z-direction one has a undisturbed diffusion and there- 
fore the mean square displacement and thus is much 
larger than in the modulated x- and y-direction. Along 
these two directions the saddles and the local maxima 
between neighboring potential valleys act as barriers for 
the dumbbell motion and therefore (equal to Dy) is 



smaller than D^. Moreover, for a dumbbell in a solvent 
the HI between the two beads comes into play which in 
general enhances the diffusion as can be seen by the shift 
between the solid and the dashed lines in Fig [21 
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FIG. 2: The mean square displacement of the CM of the 
dumbbell is shown along the x- and the j/-direction (lower 
lines) and along the ^-direction (upper lines) with HI between 
the beads (solid lines) and without (dashed lines). The pa- 
rameters are the same as in Fig. [1] 

The decay of the dumbbell diffusion as a function of 
the ratio between the modulation amplitude of the poten- 
tial and the thermal energy Vo/kBT is shown in Fig. [3] 
for one parameter set with HI (solid line) and without 
HI (dashed line) between the beads. The decay of is 
similar to the results described in Refs. d, 0] on the dif- 
fusion of point like particles. The difference between the 
two cases with and without HI is shrinking with an in- 
creasing modulation amplitude of the potential, because 
the higher diffusivity caused by HI becomes less impor- 
tant with increasing potential barriers. Accordingly, the 
ratio Dx/Dz between the diffusion along one modulated 
direction and the unmodulated direction is smaller with 
HI than for the case without HI, as shown in Fig. [3] 

In contrast to the diffusion of single point particles the 
diffusion of a dumbbell along one modulated direction 
also depends on the interplay between the two length 
scales, namely the bead distance b and the wavelength 
A of the periodic potential modulation. A typical func- 
tional dependence of the dumbbell diffusion on the ratio 
X/b is shown in Fig. |4l where the diffusion is remarkably 
enhanced for A close to Ai = 35/2. Further beyond this 
value the dumbbell diffusion decreases with increasing 
values of the wavelength up to a minimum Dx{X2), which 
is at about A2 ~ 66 for the given parameters (not shown 
in Fig. 2]). The decay of Dx in the range Ai < A < A2 can 
be explained in the following way. For increasing values 
A > Ai the beads become essentially caged within one 
single potential valley and an escape from such a trap 
gets more and more unlikely (see right inset in Fig. 3]). If 
the wavelength A is increased further beyond A2 the mean 
square displacement and the dumbbell diffusion increase 
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FIG. 3: The diffusion D^, 
ized by the free diffusion in z-direction and plotted as a 
function of Vb/ftsT for the case with HI (solid line) and with- 
out HI (dashed line). The following parameters have been 
used: fc = 10 and A = 2b. 



again and finally approach the value of the unmodulated 
case, because for large wavelengths the dumbbell does 
not feci the potential anymore. 




range the two beads hardly fit into one single potential 
valley as indicated by the middle inset in Fig. U] More- 
over, for a rather stiff spring both beads can not reach 
the minima of two neighboring potential valleys simulta- 
neously. So the required excitation energy is smaller than 
Vo and for this reason the diffusive motion is enhanced. 

In the range A < 6 the dumbbell easily finds an ori- 
entation in the potential plane where the two beads are 
located in two distinct minima, cf . the left inset of Fig. H) 
Only one bead needs to be kicked to another valley in or- 
der to make progress for the CM of the dumbbell. There- 
fore, irrespective whether the two beads belong to next 
nearest neighbor valleys or not, the required excitation 
energy for a shift of the CM depends only weakly on A. 
This is the origin of the small variations of Dj. in the 
range X < b. 




FIG. 5: The ratio between the diffusion constants D^/Dz of 
the dumbbell is shown for three different values of the spring 
constant (solid line: k = 0.1, dashed-dotted: k = 1, dashed: 
k = 10) as a function of the ratio X/b. The potential ampli- 
tude is Vo = 2fcsT. 
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FIG. 4: The normalized diffusion constant Dx/D^ {— Dy/Dz) 
of the dumbbell is shown as a function of the ratio between the 
equilibrium bead distance 6 and the wavelength A for the case 
with HI (solid line) and without HI (dashed line). The insets 
illustrate possible locations of the dumbbell with respect to 
the periodic potential for different values of X/b. The param- 
eters are Vo/fcsT = 2 and k = 10. 

If on the other hand the modulation wavelength is re- 
duced below A2, the distance from the potential mini- 
mum to a saddle is shortened. In such a case it becomes 
more likely that the dumbbell is kicked over a saddle to 
a neighboring potential valley. So the diffusion Dx in- 
creases with decreasing values in the range Ai < A < A2. 

Near the maximum of the diffusion constant at A ~ Ai 
an additional effect comes into play which supports a 
dumbbell movement from one potential valley to another 
and accordingly enhances the dumbbell diffusion. In this 
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FIG. 6: The distribution p(ri2) of the bead distances is 
shown for two different values of the spring constant, k — 10 
(solid line) and k — 1 (dashed line) and for the parameters 
Vo/ksT = 2 and X/b = 3/2. 
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The explanation given for the local maximum of the 
diffusion constant in the range of Ai in Fig.|3]is supported 
by the influence of the spring constant k on the height 
of Dx{Xi) and on the mean distance (r'12}. The local 
maximum of is especially pronounced in the case of a 
rather stiff dumbbell (see Fig. [5]) where the distribution of 
the bead distance p{r 12) is not changed by the periodic 
potential (see Fig. On the other hand for smaller 
values of k the maximum of /o(r 12) is more and more 
shifted from & to A. In this case the beads relax down 
to the potential valleys, so a higher excitation energy is 
required for a shift of the CM, which results in a smaller 
diffusion constant as indicated by the solid line in Fig. [5] 
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FIG. 7: The orientational distribution p(4') of the dumbbell 
axis in the xy-plane is shown for two different values of the 
spring constant (solid line: = 1, dashed: k = 10) and for 
Vo/fesT = 2. The angle <E> is measured with respect to the 
a;-axis and the ratio between the wavelength A and the bead 
distance h is in part a) \/b = 2 and in part b) X/b — 1. 

The distance between the potential valleys depends on 
the direction in the xy-planc and therefore the orienta- 
tional distribution of the dumbbell axis p{^) in Fig. [3 
provides a complementary information to p{ri2). The 



two beads of the dumbbell may relax more easily to the 
potential minima in the case of a soft spring (cf. solid 
lines in Fig. [7]) compared to a stiff spring (cf. dashed 
lines). So the orientational distribution of the dumbbell 
axis in the a;?;-plane shows, besides the maxima along 
the X- and the y-direction, a local maximum for a diag- 
onal orientation of the dumbbell axis. This is displayed 
in part a) of Fig. [71 If the wavelength is reduced to 
X — b one finds, in addition to the maximum in the di- 
agonal direction (11), local maxima along the (21)- and 
the (12)-direction, as indicated by part b) of Fig. [71 

Conclusions. - The Brownian motion of a dumbbell in a 
two-dimensional periodic potential has been investigated 
in terms of a Langevin model. For an increasing barrier 
height between the potential minima, we find a decreas- 
ing diffusion constant for the center of mass along the 
modulated spatial directions as well as a reduction of 
the influence of the hydrodynamic interaction. For stiff 
springs the diffusion additionally depends on the ratio 
between the wavelength A of the potential modulation 
and the equilibrium dumbbell extension b. In the range 
A ~ 36/2 this interplay is especially pronounced, because 
the two beads do not fit into one single or two neighboring 
potential minima anymore and this mismatch causes a re- 
duced effective barrier height and thus an enhanced dif- 
fusion constant. If the spring constant is small the beads 
can relax down to the potential minima over a wide range 
of A, which results in a diffusion constant that is rather 
independent from A in this domain. So the height of the 
maximum of the diffusion constant at A w 3&/2 increases 
with the spring stiffness. For modulation wavelengths 
further beyond 36/2 the diffusion constant decays mono- 
tonically until some minimum is reached. In this range 
the dumbbell is essentially caged in one single potential 
valley and it is rather unlikely that it escapes. Beyond 
this minimum as A goes to infinity the diffusion constant 
grows until the free diffusion limit is reached. 

The presented results may be useful for sorting poly- 
disperse particle mixtures with respect to the particles' 
elasticity and size. 
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